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A new method called difFusion factorial moment (DFM) is used to obtain scaling features em- 
bedded in spectra of complex networks. For an Erdos-Renyi network with connecting probability 
Per < jf, the scaling parameter is 5 = 0.51, while for per > jj the scaling parameter deviates 
from it significantly. For WS small- world networks, in the special region pr £ [0.05,0.2], typical 
scale invariance is found. For GRN networks, in the range of S £ [0.33, 049], we have 5 = 0.6 ± 0.1. 
And the value of 5 oscillates around 5 — 0.6 abruptly. In the range of 6 £ [0.54, 1], we have basically 
5 > 0.7. Scale invariance is one of the common features of the three kinds of networks, which can 
be employed as a global measurement of complex networks in a unified way. 
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I. INTRODUCTION 

In recent years, complex networks attract special at- 
tentions from diverse fields of research [1]. Though 
several novel measurements, such as degree distribu- 
tion, shortest connecting paths and clustering coeffi- 
cients, have been used to characterize complex networks, 
we are still far from complete understanding of all pe- 
culiarities of their topological structures. Finding new 
characteristics is still an essential role at present time. 

Describe the structure of a complex network with the 
associated adjacency matrix. Map this complex network 
with N nodes to a large molecule, the nodes as atoms and 
the edges as couplings between the atoms. Denote the 
states and the corresponding site energies of the atoms 
with , |2) , ... jA^)} and {ei,e2, ■ ■ ■ En}, repsectively. 
Consider a simple condition where the Hamiltonian of 
the molecule reads, 

Hfrin — (-^o)mn ~t~ {^QQupi^ng^mn: (1) 

where, 

mn : 

,„ N _ J 1 — 5mn connected (2) 

{tlcoupUng)mn - | q disconnected 

By this way a complex network is mapped to a quan- 
tum system, and the corresponding associated adjacency 



matrix to the Hamiltonian of this quantum system. 

The structure of a complex network determines its 
spectrum. The characteristics of this spectrum can re- 
veal the structure symmetries, which can be employed as 
global measurements of the corresponding complex net- 
work [2-12]. In our recent papers [13-16], several tempo- 
ral series analysis methods are used to extract character- 
istic features embedded in spectra of complex networks. 

In the present paper, a new concept, called diffusion 
factorial moment (DFM), is proposed to obtain scale fea- 
tures in spectra of complex networks. It is found that 
these spectra display scale invariance, which can be em- 
ployed as a global measurement of complex networks in 
a unified way. It may also be helpful for us to construct 
a unified model of complex networks. 

II. DIFFUSION FACTORIAL MOMENT (DFM) 

Represent a complex network with its adjacency ma- 
trix: A{G). The main algebraic tool that we will use for 
the analysis of complex networks will be the spectrum, 
i.e., the set of eigenvalues of the complex network's ad- 
jacency matrix, called the spectrum of the complex net- 
work, denoted with {Em \m — 1, A^ }. Connecting the be- 
ginning and the end of this spectrum, we can obtain a 
set of delay register vectors as [13], 

{El — Eq,E2 — El, ... En — En-l} 
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{E2 — Ei,E^ — i?2, . . . En+l — En} 



2 



{E 



N 



E 



JV-l 



N, ■ • 



■ En-2 — -Era-s} 



{Eq — En, El — Eq, . . . En-i — En-2} 



(3) 



Considering each vector as a trajectory of a particle 
during n time units, all the above vectors can be regarded 
as a diffusion process for a system with N + 1 particles 
[17]. Accordingly, for each time denoted with n we can 
reckon the distribution of the displacements of all the 
particles as the state of the system at time n. 

Dividing the possible range of displacements into 
Mq bins, the probability distribution function (PDF) 
can be approximated with ~ , where 

Km{n) \m = 1,2, ...Mq is the number of particles whose 

displacements fall in the m'th bin at time n. To obtain 
a suitable Mq, the size of a bin is chosen to be a fraction 

of the variance, s - 



^k=l (^k—Ek^l)^ 



N+l 

If the series constructed with the nearest neighbor level 
spacings, {i^i — Eq, E2 — Ei, . . . Eq — En}, is a set of ho- 
mogeneous random values without correlations with each 
other, the PDF should tend to be a Gaussian form when 
the time n becomes large enough. Deviations of the PDF 
from the Gaussian form reflect the correlations in the 
time scries. Here, we are specially interested in the scale 
features in spectra of complex networks. 

Generally, the scale features in spectra of complex net- 
works can be described with the concept of probability 
moment (PM) defined as [18], 



Mn 



(4) 



m=l 



where Pm is the probability for a particle occurring in the 
m'th bin. Assume the PDF takes the form, 



Pm{n) = X-F (^) ■ 



An easy algebra leads to. 



In C,(n) = A + 5-{l-q) ln(n) 



(5) 



(6) 



If the considered series is completely uncorrelated, the 
resulting diffusion process will be very close to the con- 
dition of ordinary diffusion, where 5 = 0.5 and the func- 
tion (^) in the PDF is a Gaussian function of 
5 ^ 0.5 can reflect the departure of the diffusion process 
from this ordinary diffusion condition [19]. The extreme 



condition is the ballistic diffusions , whose PDFs read 
Pm{n) = F {^)- The values of 5 at this condition are 
1. 



But the approximation of PDF, p„ 



the 



above computational procedure will induce statistical 
fluctuations due to the finite number of particles, which 
may become a fatal problem when we deal with the spec- 
trum of a complex network. The dynamical information 
may be merged by the strong statistical fluctuations com- 
pletely. Capturing the dynamical information from a fi- 
nite number of cases is a non-trivial task. 

This problem is firstly considered by A. Bialas and R. 
Peschanski in analyzing the process of high energy colli- 
sions, where only a small number of cases can be avail- 
able. A concept called factorial moment (FM) is pro- 
posed to find the intermittency (self-similar) structures 
embedded in the PDF of states [18,20-24]. The definition 
of FM reads. 



F,{M) 



M 
m=l 



{Jm-l)---{Jm-q + ^), 



(7) 



where M is the number of the bins the displacement 
range being divided into and Jm the number of parti- 
cles whose displacements fall in the m'th bin. 

Stimulated by the concept of FM, we propose in this 
paper a new concept called diffusion factorial moment 
(DFM), which reads. 



DFMg{n) 



Mo 



{K^-l)---{Krn-q+l). (8) 



m=l 



Herein we present a simple argument for the ability of 
DFM to filter out the statistical fluctuations due to fi- 
nite number of cases [20, 21]. The statistical fluctuations 
will obey Bernoulli and Poisson distributions for a system 
containing uncertain and certain total number of parti- 
cles, respectively. For a system containing uncertain total 
number of particles, the distribution of particles in the 
bins can be expressed as. 



Q {Ki, K2,--- , Kmo \pi,P2--- ,PMo ) 

~ K^\K2\-- Kmo^-^'^ ^2 PMo ' 

where K = Ki+ K2 + h Kmo- Hence, 

{K,n{K,n-l)---{K^~q+l)) 
= J dpi_dp2 ■ ■ ■ dpM„P{Pl,P2, ■ ■ ■ ,PMo)- 
EE-- - E Q{Ki,K2,--- ,Kmo \Pl,P2, 



(9) 



Ki K2 



,PMo) 



Mo 



KmiKra - I) ■ ■ ■ {Krr, - q + I) 

= K{K - 1){K -2)---{K-q + l)J dpidp2 ■ ■ ■ dpuo' 

P{P\,P2,--- ,PMo)-PIi 

= K{K-l)...{K-q+l){p<Q 

(10) 
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That is to say, 

DFMq{n) oc Cq{n) 
And consequently, Eq.{6) becomes, 



(11) 



InDFMqin) = C + lnC,(n) = B + S ■ {1 - q)lii{n) (12) 

Therefore, DFM can reveal the strong dynamical fluctua- 
tions embedded in a time series and filter out the statisti- 
cal fluctuations effectively. We will use the DFM instead 
of the PM to obtain the scale features in spectrum of a 
complex network. 

It should be pointed out that the scale features in our 
DFM is completely different from that in FM. The FM re- 
veals the self-similar structures with respect to the num- 
ber of the bins the possible range of the displacements 
being divided into, i.e., the scale is the displacement. In 
DFM, the considered scale is the time n. At time n, 
the state of the system is, (i?„ — Eq, E^+i — Ei,En+2 — 

E2, • • • , -Eat — Ejy^q, Ef) — En-q+i)- 

In one of our recent works [13], joint use of the dc- 
trended fluctuation approach (DFA) and the diffusion 
entropy (DE) is employed to find the correlation features 
embedded in spectra of complex networks. In that paper 
we review briefly the relation between the scale invariance 
exponent, S, and the long-range correlation exponent a. 
For fractional Brownian motions (FBM) and Levy walk 
processes, wc have 6 = a and S — ^^2-0 ' respectively. 
Generally, we can not derive a relation between these two 
exponents. Herein, we present the relation between the 
concepts of DFM and DE. Prom the probability moment 
in EqA we can reach the corresponding Tsallis entropy, 
Srsaiiis, which reads. 



JTsaUis 



Ma 



1-q 



1-q 



(13) 



A trivial computation leads to the relation between the 
DE, (denoted with Sde), the PM and the Tsallis entropy, 
as follows, 



S 



DE 



lim - 
limi^ 



Mo 

E 



= lini STsallisiq)- 



(14) 



Hence DFM can detect multi-fractal features in spectra 
of complex networks by adjusting the value of q. The DE 
is just a special condition of DFM with g — » 1. What is 
more, the DFM can filter out the statistical fiuctuations 
due to finite number of eigenvalues in the spectrum of a 
network. 

The adjacency matrices arc diagonalized with the Mat- 
lab version of the software package PRO PACK [25] . 
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FIG. 1: Four typical results for the Erdos-Renyi network 
model. The parameter q = 3. Denote the size of a network 
with N. {a.)N = W*, per = ^ < Pc- We have 5 = 0.51, 
which is consistent with the random i)ohavior of the spectrum. 
Tiio corresponding PDF is Gaussian, (b) N = 10*, Per = 
— pc- Wc have S = 0.60, a slight deviation from random. 
For (c) and (d) {N,per,S) = (4 x 10^, 0.68) and (4 x 
10^, |-,0.87), respectively. 



III. RESULTS 

Consider firstly the Erdos-Renyi model [26,27]. Start- 
ing with N nodes and no edges, connect each pair with 
probability Per- For per < ^/n the network is broken 
into many small clusters, while ior Per > Vjv a large clus- 
ter can be formed, which in the asymptotic limit contains 
all nodes [27]. Pc — ^/n is a critical point for this kind of 
random networks. 

Fig.l presents four typical results for Erdos-Renyi net- 
works. For Per < Vjv, the scaling exponent is, (5 = 0.51, 
which is consistent with the random behavior of the spec- 
trum. With the increase oi per, S becomes larger and 
larger. The spectrum tends to display a significant scale 
invariance. 

As one of the most widely accepted models to cap- 
ture the clustering effects in real world networks, the WS 
small world model has been investigated in detail [1,28- 
31]. Here we adopt the one-dimensional lattice model. 
Take a one-dimensional lattice of L nodes with periodic 
boundary conditions, and join each node with its k right- 
handed nearest neighbors. Going through each edge in 
turn and with probability Pr rewiring one end of this edge 
to a new node chosen randomly. During the rewiring pro- 
cedure double edges and self-edges are forbidden. 

Fig. 2 and Fig. 3 show the results for two extreme con- 
ditions of the WS network model, i.e., the regular net- 
works with different right-handed neighbors {pr = 0) and 
the corresponding completely rewired networks (pr = 
1). When the value of k is unreasonable large (fc = 
20, 40, 80, 120), the DFM wifl not obey a power-law. The 
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FIG. 2: DFM for the two extreme conditions of the WS 
network model, i.e., regular networks (pr ~ 0) and the corre- 
sponding completely random networks {pr = 1). The size of 
a network N = 3000. And q — 3. For these generated net- 
works, when the number of right-handed neighbors k is small 
(fc = 1,2,- • -,9) the DFMs obey a power- law. 



FIG. 4: Four typical results for WS small-world model. The 
parameters q = 3, fe = 2. The size of a network is 3000 . 
(a) Pr = 0.0, S = 0.67; (b) pr = 0.1, S = 0.71; (c) pr = 0.3, 
S = 0.64; {d)pr = 0.8, S = 0.80. 
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FIG. 3: The values of the exponent 5 for the two extreme 
conditions of the WS network model, i.c, the regular networks 
{pr ~ 0) and the corresponding completely random networks 
{pr = 1). The size of a network N = 3000. And g = 3. 
The values of S for the regular networks arc in the range of 
0.58 ± 0.08, a slight deviation from 0.5 that corresponds to 
a Gaussian distribution. The values of 5 for the completely 
random networks are significantly larger than that of the cor- 
responding regular networks (with a few exceptions). 



scaling exponents for the regular networks are basically 
in the range of ^ = 0.58 ± 0.08, a slight deviation from 
that of the Gaussion distribution. The scaling exponents 
for the completely rewired networks with k = 2, 3,4, 6, 7 
are significantly larger than that of the corresponding 
regular networks. 

Four typical results for the networks generated with 



FIG. 5: The values of 5 for generated WS networks with 
different rewiring probabilities. The parameters q = 3, k = 2. 
In the special range of pr £ [0.05,0.2], where the WS small 
world network model can capture the chaxacteristics of real 
world complex networks, we have 5 £ 0.71 ± 0.05. 



the WS model with different rewiring probability values, 
as shown in Fig. 4, illustrate the significant scale invari- 
ance in spectra of these WS networks. The values of 6 for 
these generated networks with k = 2 and k = 5 are pre- 
sented in Fig5 and Fig. 6, respectively. We arc specially 
interested in the rough range of Pr G [0.05, 2] where the 
WS model can capture the characteristics of real world 
networks. For the generated networks with fc = 2, in the 
range of pr G [0.05, 0.2] we have 6 G 0.71 ± 0.05. And in 
the condition of fc = 5, ^ is 0.85 ± 0.05 in the range of 
Pr e [0.1,0.2]. 

Consider thirdly the growing random network (GRN) 
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FIG. 6: The values of 5 for generated WS networks with 
different rewiring probabiUties. The parameters g = 3, fc = 5. 
In the special range of pr € [0.1,0.2], where the WS small 
world network model can capture the characteristics of real 
world complex networks, we have S € 0.85 ± 0.05. 
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FIG. 7: Four typical results for networks generated with the 
GRN model. The parameter q = 3. The size of a network 
is 4000 . {a) e = 0.0, S = 0.88; (b) 9 = 0.33, S = 0.53; (c) 
e = 0.5, d = 0.68; {d)9 = 1.0, d = 0.86. 



model [30,32]. Take several connected nodes as a seed. 
At each time step, a new node is added and a link to one 
of the earlier nodes is created. The connection kernel Ah, 
defined as the probability that a newly introduced node 
links to a pre-existing node with k links, determines the 
structure of this network. The considered complex net- 
works are generated with a class of homogeneous connec- 
tion kernels, ^fe oc fc^(0 < 61 < 1). 

The arguments in literature [32] show that there are 
two critical points at = ^ and 62 = ^, which sepa- 
rate the networks into four groups. The four groups are 
[0, 5),5and (i, 1). From the values of 6 for GRN 
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FIG. 8: The values of 5 for the networks generated with the 
GRN model with different 6. At two points 6 = 0.33, 0.49, we 
have S — 0.53 and 0.52 (two minimum values), respectively. 
In the range of 61 £ [0.33, 0.49], we have 5 = 0.6 ± 0.1. And 
the value of 6 oscillates around 5 = 0.6 abruptly. In the range 
of € [0.54, 1], we have basically 5 > 0.7. 



networks with different 6, shown in Fig. 6, we can find 
that at two points 6 ~ 0.33, 0.49, wo have S = 0.53 and 
0.52 (two minimum values), respectively. In the range of 
e e [0.33,0.49], we have 5 = 0.6 ± 0.1. And the value 
of S oscillates around S = 0.6 abruptly. In the range of 
G [0.54, 1], we have basically d > 0.7. 



IV. SUMMARY 

In summary, we introduced a new concept called DFM 
and use it to reveal scale invariance features embedded 
in spectra of complex networks. For an Erdos-Renyi net- 
work with connecting probability psn < j^- the scaling 
exponent is S = 0.5, while for per > jj the scaling ex- 
ponent deviates from 0.5 significantly. For the regular 
networks generated with the WS model with pr = 0, 
the scaling exponents deviate slight from 0.5, the value 
corresponding to the Gaussian PDF. The other extreme 
condition is that the S values for the random networks 
generated with the WS model with pr ~ I are basically 
significant larger than that for the corresponding regu- 
lar networks (there are few exceptions). In the specially 
interested range of G [0.05,0.2], where the WS model 
can capture the properties of real world networks, the 
spectra display a typical scale invariance. Two critical 
points are found for GRN (growing random network) 
networks at = 0.33 and 9 ~ 0.49, at which we have 
two minimum values of 6 = 0.53,0.52, respectively. In 
the range of ^ e [0.54,1], we have basically 6 > 0.7. 
Hence wc find self-similar structures in all the spectra 
of the considered three complex network models. This 
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common feature may be used as a new measurement of 
complex networks in a unified way. Comparison with 
the regular networks and the Erdos-Renyi networks with 
Peb^ < Pc — jj tells us that this self-similarity is non- 
trivial. 

The self-similar structures in spectra shed light on 
the scale symmetries embedded in the topological struc- 
tures of complex networks, which can be used to obtain 
the possible generating mechanism of complex networks. 
Quasicrystal theory tells us that the aperiodic structure 
of lattice will induce a fractal structure in the correspond- 
ing spectrum. The most possible candidate feature shar- 
ing by all the complex networks constructed with the 
three models may be fractal characteristic, which has 



been proved in a very recent paper [33]. Based upon 
this feature, we may construct a unified model of com- 
plex networks. 
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